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The first locations of short gamma-ray bursts (GRBs) in elliptical galaxies 
suggest they are produced by the mergers of double neutron star (DNS) bina¬ 
ries in old stellar populations. Globular clusters, where the extreme densities 
of very old stars in cluster cores create and exchange compact binaries effi¬ 
ciently, are a natural environment to produce merging NSs. They also allow 
some short GRBs to be offset from their host galaxies, as opposed to DNS 
systems formed from massive binary stars which appear to remain in galactic 
disks. Starting with a simple scaling from the first DNS observed in a galactic 
globular, which will produce a short GRB in ~300My, we present numeri¬ 
cal simulations which show that ~10—30% of short GRBs may be produced 
in globular clusters vs. the much more numerous DNS mergers and short 
GRBs predicted for galactic disks. Reconciling the rates suggests the disk 
short GRBs are more beamed, perhaps by both the increased merger angular 
momentum from the DNS spin-orbit alignment (random for the DNS systems 
in globulars) and a larger magnetic field on the secondary NS. 

This paper has been submitted to Nature Physics. In accordance with the editorial policy 
of Nature, it is embargoed for discussion with the press. Any questions should be directed 
to Josh Grindlay, jgrindlay@cfa.harvard.edu. 

The origin of cosmic gamma-ray bursts, like most astronomical mysteries, required 
precise source positions and identification at other wavelengths for the first breakthroughs. 
After nearly 25 years, the so-called “long” GRBs (^2-200sec), were located by coded 
aperture imaging of their hard X-ray emission which enabled more precise positions from 
their soft X-ray afterglows^ and then optical counterparts^ to be found. The optical 
identifications established them to be at cosmological distances. They are now understood 
to be due to relativistic jets produced in the core collapse of massive stars to form stellar 
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mass black holes in certain ’’hyper” supernova events.^ However, about a third of GRBs 
were previously recognized"^ to be distinctly different: both shorter duration (^0.2-2 sec) 
and harder spectra. The recently-launched Swift satellite® was designed to solve the 
origin of the short GRBs by rapid detection of their X-ray and/or optical afterglows to 
enable their identihcation. 

The hrst of these short GRBs, GRB050509b, was detected and located precisely (^5 ” ) 
enough by Swift^ to enable its plausible optical identification^ with the halo of an elliptical 
galaxy at 1.12 Gpc distance. A second short burst, GRB0507024, was located precisely 
by Swift^’^ with an optical counterpart detected at an offset of ~2.5kpc from the center 
of another elliptical (E2) galaxy at very similar redshift (z = 0.257) and hence distance. 
Since elliptical galaxies have long ceased active star formation, this suggested^^ that short 
GRBs are associated with old stellar populations and likely due to the merger of two 
neutron stars or a neutron star (NS) and stellar black hole (BH) in a compact binary 
system as originally suggested for GRBs generally.The offset position for GRB050509b 
at 40±13kpc from its elliptical galaxy Gl^ would previously be attributed^^ to the ejection 
of the progenitor compact binary from the galaxy by the kick(s) imparted to it by the 
supernova event(s) which created its constituent neutron star(s) or black hole. This 
implies a kick velocity ^600 km s“^ to remain bound to the elliptical, and yet ^200 km 
s“^ to have its apogalacticon in the halo. However, this is inconsistent with it being like the 
Hulse-Taylor binary pulsaror the six other DNS systems now known in our Galaxy that 
were produced from the evolution of a massive binary system^"^’^® since recent analysis^® 
shows these are all likely low-kick velocity (^50 km s“^) systems and thus are expected to 
remain within the central potential of their parent galaxies. Thus short GRBs from spiral 
or star formation galaxies, where massive binary evolution is producing DNS systems, are 
expected in their disks. Since the third short GRB located by Swift^'^ is also associated 
with an elliptical, an alternative origin for their DNS progenitors is suggested. 

We suggest instead that short GRBs are at least partly produced by DNS systems 
which formed in globular cores by exchange interactions of NSs into previously formed 
compact binaries containing a NS and low mass secondary. Subsequent mergers of these 
dynamically formed DNS compact binaries produced the short GRBs located by Swift. 
There is direct evidence that DNS systems form by dynamical interactions in globular 
clusters: the millisecond pulsar M15-G with a NS binary companion in the core collapsed 
globular cluster M15 in our Galaxy^® must have formed within its ~100 My spindown 
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age and then was scattered out of the cluster core by a NS exchange interaction^® with a 
cluster low mass X-ray binary (LMXB) to produce a DNS with present eccentricity e = 
0.68 and binary period P = 0.33d. We demonstrate this formation scenario, and hnd that 
signihcantly more DNS systems are produced by NSs in core collapsed globular clusters 
(like M15) by the exchange of a NS for the lower mass secondaries in quiescent low mass 
X-ray binary (qLMXB) or millisecond pulsar (MSP) binary systems in the cluster cores. 
A previous scenario^® for long GRBs from globular clusters, brought to our attention 
after hnal submission of this paper, invoked NSs “smothered” in stellar merger collisions 
to form BHs (for which cross sections remain unknown) but produced comparable rates 
to our detailed results for re-exchange production of DNS systems. 

Recent Chandra and HST observations^^’^^ have shown that MSPs in globular cluster 
cores do indeed undergo subsequent encounters to re-exchange their evolved low mass 
secondaries for a main sequence (MS) star, which is not found or expected^"^ for isolated 
evolution from a LMXB or qLMXB. The 30ms spin period of the M15-C pulsar implies it 
did not complete its spinup (90% of MSPs in globulars have shorter periods^^), so that its 
original companion, as a LMXB, was thus a ~O.4-O.7M0 MS star (most probable in the 
core) that was exchanged for a cluster NS, in many cases another MSP. The re-exchange 
formation model developed here can operate with a wide range of original secondaries 
(most often ~0.05-0.2Mq He white dwarfs). 

Given the existence of M15-G, can the globular cluster populations of elliptical galaxies 
plausibly contain enough DNS systems to account for a signihcant fraction of the short 
GRB rate? The answer is yes. We start with the conservative assumption that each 
L* galaxy contains a globular cluster system with 10^ globulars (typical for L* ellipticals) 
and that the number density of such galaxies is ~0.01 Mpc“^ as determined^^ from 2MASS 
and 2dF galaxy counts, yielding a globular cluster space density ncc lOMpc Then, 
using the “local rate” of observed short GRBs derived by Schmidt^® as Rg_rb ~0.08 Gpc“^ 
yr“^, the “observed” fractional number per globular of DNS systems like M15-G, with its 
Tmerge ~300My merger time, would need be Gfes = RcRsTmergel'^GC ~2.4 x 10“^ to account 
for short GRBs. Gorrecting for the short GRB beaming factor of 30-50 inferred^® for 
(only!) one short GRB, GRB050724, multiplies the observed value for the total number 
of short GRB progenitors. However this is partly offset by the pulsar beaming factor for 
MSPs,^® Bmsp %2, where the lower limit applies if DNS systems contain at least one MSP. 
Thus the net beaming factor Bdns = Bgrb/Bmsp ^15 and the total number of DNS 
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systems needed per globular is foNS = fobs ^dns ^0.04. Thus short GRBs could then 
arise from a population of DNS systems in globular clusters around L* galaxies if they 
were produced within ~4% of its globulars. If, as we show below, this largely arises in the 
~20% of globulars that have undergone core collapse,then DNS systems are required 
in ~20% of the post core collapse (PCC) clusters, which would suggest another ~7 of the 
~40 PCC globulars in the Calaxy should contain a M15-C-hke system. This is consistent 
with currently incomplete globular cluster surveys for pulsars which already hnd 8 MSPs 
in binaries with large eccentricities,^^ many of which could have NS companions. 

Motivated by this simple scaling, we have conducted numerical scattering experiments 
using the scatters and sigmaS programs^® in the Starlab software environment^® to 
compute cross sections for the formation of DNS systems with merger times less than 
the age of the universe, as well as the distributions of their emergent velocities, retention 
probabilities, hnal periods and eccentricities. Each experiment consists of a three-body 
scattering between a target binary containing a 1 . 4 M 0 NS primary and a MS or helium 
WD secondary, with parameters chosen to represent the MSP, qLMXB and (relatively 
rare) LMXB targets that a lone NS in the cluster might encounter. Our target outcome 
is an exchange interaction leading to the formation of a DNS binary with a merger time 
smaller than a typical cluster age, ~ 10 Cyr. 

We adopt the following parameters for our experiments: MS secondaries are assumed 
to have masses of 0.2, 0.4, or O.SMq, while the WD masses are 0.05 or O. 2 M 0 . The 
initial binary orbits have periods ranging from 0.01 to lOOd (as appropriate for WD or 
MS secondaries) and zero eccentricity, as observed (very nearly) for most MSPs (and 
as expected for their progenitor qLMXBs). Incoming NS velocities are 10 km s“^, a 
typical cluster velocity dispersion. We include the possibility of stellar collisions during 
the interaction. The secondary radii are the zero-age MS radii for the MS stars, and 0.03 
and 0.018 i ?0 for the 0.05 and O. 2 M 0 WD secondaries, respectively. NS radii are taken to 
be 15 km. We performed a total of ~ 3 x 10® simulations, leading to statistical errors in 
our derived cross sections of .<3%. 

The results are quite remarkable: M15-C systems could be produced with cross sections 
that are comparable (~0.8 AU®) over target binary periods Pb ~0.1-ld (Figure 1). Even 
more striking is the range of output DNS periods (Pb) and eccentricities (e) produced by 
these exchange encounters (Figure 2). For the O. 4 M 0 MS secondary case appropriate 
to M15-C, and the exchange interaction randomly distributed between t = 0 and t = 
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10^° years, the orbital evolution due to the emission of gravitational wave radiation is 
subsequently followed for each of the DNS systems produced until the current age of the 
Galaxy of 10^'^ years. The star indicates the observed position of the binary pulsar M15-C 
and is “predicted” by the distribution of survivors (heavy dots). 

With the cross sections in Fig 1 we can estimate the number of DNS systems produced 
per globular. We assume that the orbital period distribution of progenitor systems is 
similar to that observed for MSPs in Galactic globular clusters,resulting in an average 
cross section of a ~ 0.8 AU^. We note that this greatly exceeds the area of the target 
binary orbit due to gravitational focusing at low encounter velocity ( r'U 10 km s ^). The 
instantaneous formation rate of DNS binaries which will merge within 10 Gyr is then given 
by T = Npj-nav, where Apr is the number of potential progenitor systems in a cluster with 
mean neutron star number density n and velocity dispersion v. After substitution and 
scaling to parameters typical of globular clusters, we obtain 

^ (p^) (lol^) (llok^) (m) 

where, as discussed below, we normalize to a neutron star density of n ~ 10® pc“^ and 
a total of Apr = 20 progenitor binaries containing a NS primary. The latter hgure 
is estimated from the population of target LMXBs (~0.05), qLMXBs (~1) and MSPs 
(~10) per cluster, as derived from Chandra studies of (primarily) non-core-collapse but 
high central density clusters,^® but increased by a factor of 2 for the expected production 
during core collapse of additional target LMXBs, MSP-WD and MSP-MS systems. 

The quantity Tdns varies dramatically over a cluster’s evolutionary lifetime, rising 
sharply around the time of core collapse. Assuming a 10% neutron-star retention frac¬ 
tion,®^ we estimate that the Galaxy’s “non-core-collapse” clusters have n ~ 10^ —10^pc“®, 
while models for core-collapse clusters®^ predict n ^ 10® pc“®. The cluster central velocity 
dispersion v depends only weakly on the state of the core,®^ scaling as n®'®®. As a result, 
during the homologous late stages of core collapse the central density n scales roughly as 
where T = tcc — tis the time remaining until collapse, and the total number of DNS 
systems produced in the cluster, Ndns = /hDArshf, is dominated by the behavior near 
r = 0. Thus post-core-collapse (PGG) clusters, ~20% of the total in the Galaxy,^^ would 
dominate the overall production of DNS binaries in globulars. 

We use the M15 simulation data of Dull et al.®^ to resolve the formal singularity at 
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T = 0, and find that 


n(t) dt ^ 20noro , 


( 2 ) 


^ ^cc — 

where Tq is the time remaining until core collapse from a state with central density 
no = n{tcc — To). The dependence of n on r makes the precise instant at which this 
expression is evaluated relatively unimportant. We choose no = 10®pc“^, corresponding 
to To lOOtfc ~ lOOMyr for parameters appropriate to M15. The behavior of clusters 
during the PCC phase remains poorly understood but simulations indicate that the re¬ 
expansion is signihcantly slower than the collapse, and may include a series of nonlinear 
core oscillations. Thus we expect that the post-collapse phase will contribute at least 
as much to Ndns as did the collapse; each recollapse will contribute roughly as much 
again. Weighting the cross sections of the DNS binaries with the observed orbital period 
distribution for MSPs in globulars,^^ since they dominate the DNS formation, we then 
estimate the total number of merging DNS binaries per PCC cluster as 


A^dns ^12 
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The coefficient 12 in Eq. (3) comes from ~ 25% of the progenitors scattering directly 
into merging orbits, and a further ~ 35% forming wider DNS systems that harden into 
merging orbits following a subsequent interaction during the core-collapse phase. We note 
that ~10% of these systems will stem from (q)LMXBs and hence should be comparable 
to M15-C. The remainder, descended from MSPs, will be wider and more eccentric. 

Thus, combining the two populations (qLMXBs/LMXBs and MSPs), ignoring the 
small differences in donor mass in qLMXBs and LMXBs, we derive a total of ~12 DNS 
binaries per core-collapse cluster, or ~ 480 DNS systems in the Galactic system of ~40 
PCC globulars formed over the age of the Galaxy. Some 400 of these have already merged 
due to the emission of gravitational waves (we estimate that in the last 1 Gyr ~40 systems 
have merged), and ~80 survive today (but will merge within a Hubble time, T/j). The 
majority of DNS systems are expected to remain in the globular cluster, as the exchange 
encounters tend to induce only a small recoil velocity on the binary pulsar. Only the 
LMXBs have signihcant (~40%) ejection of DNS pulsars into the Galactic disk and halo, 
as suggested for M15-C.^® Our derived total of ~80 DNS systems surviving in the globular 
cluster system of the Galaxy (i.e., ~2 per PCC globular) that will merge within T/^ and 
produce short GRBs is ~ lOx larger than our simple scaling rate from M15-C alone. The 
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difference comes from the fact that M15-C itself has a shorter than average merger time; 
our calculations now sample the entire distribution of merger times. 

The derived DNS merger rate in the Galaxy of 40 per Gyr in 200 (total) globular 
clusters implies, for the same scalings to globular clusters in external galaxies, ~2 DNS 
mergers per year per Gpc^. With the estimated short GRB beaming factor^*^ Bgrb, 
this gives an observable GRB rate of ~0.07 Gpc“^ or roughly the Schmidt rate^® 
but a factor of ~1.5-100 below the rates (with large uncertainties) recently obtained by 
Guetta and Piran,^^ for DNS production tied to the star formation rate. Given the large 
uncertainties in short GRB beaming factors, short GRB rates, and the details of core 
collapse (e.g., core collapse oscillations could increase our rates by factors ^2), and the 
additional DNS production from cluster primordial binaries containing NSs (most NSs 
retained in the cluster are likely in binaries), we conclude that perhaps ~10-30% of short 
GRBs are from DNS mergers in globulars. 

Our DNS merger rate in globulars in the Galaxy of ~400 in 10^° years or ~0.04Myr“^ 
yields a detection rate for ground based gravitational wave antennas that is ^200 x smaller 
than the rates summarized^®’^^ for DNS mergers in high mass binaries in the Galactic 
disk. Either these high mass progenitor rates for DNS production in the disk have been 
signihcantly overestimated, or most do not produce short GRBs, or their beaming angles 
are much smaller. We note that the NS population in globulars is likely dominated by 
recycled pulsars with low magnetic helds (B ~ 10®“® Gauss) whereas in the disk DNS 
systems the second-born NS, with no recycling, will have B rsj 10 i 2 14 Although 

this B held is inconsiquential for the merger, it may launch a magnetically driven jet 
with smaller beaming angles. Gombined with the randomly aligned spin-orbital angular 
momenta of DNS systems in globulars vs. the aligned spin-orbit vectors for disk systems, 
disk short GRBs might have beaming angles ^10 x smaller to reconcile the rates. 

The recent discovery of an interstellar medium (ISM) with density njsM ~ 0.1 cm“® 
in the globular cluster 47 Tucanae®^ from the variable dispersion measure of its radio- 
detected MSPs is consistent with the ISM density deduced^® for GRB050724. Globular 
clusters that have just traversed their host galaxy disks will have lower ISM densities, so 
that fainter afterglows are expected for some short GRBs, such as appears to be the case 
(given the non-detection) for GRB050509b.® 

Finally, the short GRB050709 discovered with HETE®^ was followed by a faint, soft 
and delayed (~150s) X-ray flare which may point to a different origin. Mergers of BH-NS 
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binaries may disrupt the NS and produce a “delayed” short GRB, whereas the lack of 
direct evidence for stellar mass BHs in globulars restricts the sample to NS-NS mergers. 
Indeed GRB050709 is reported^® to be associated with a star formation galaxy and so is 
more plausibly a DNS merger from a massive binary although of course its host galaxy 
also presumably includes globular clusters. Thus short GRBs plausibly arise from two 
populations, as suggested on other grounds^^: DNS or NS-BH mergers from high mass 
systems in star-forming disks with smaller beaming and DNS mergers from dynamically 
formed systems in globulars. The latter may be conhrmed by further identihcations with 
ellipticals or locations in host galaxy halos; the globular clusters themselves are too faint 
for detection with magnitudes V ^29 even at z ~0.1. 
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Figure 1. Cross sections for DNS production by exchange of low mass secon¬ 
daries (in MSPs or LMXBs) with mass indicated for a NS. The O.OSMq and 
O. 2 M 0 secondaries are calculated for He-WD mass-radius models and shown as the heavy 
curves, whereas the remaining ^0.2 Mq cases are for main sequence stars. Target binary 
periods span the range shown. Corresponding compact binaries represented as actual 
“targets” are: LMXBs, with MS secondaries and masses ^0.2 M©; qLMXBs, with (typ¬ 
ically) O.OSMq or 0.2MqWD or 0.2Mq main sequence secondaries; and MSPs, with 
O.OSMq or 0.2Mq WD secondaries. The cross sections are comparable ( rs^O. 8 AU2) 
over the distribution of target orbital periods (mostly ~0.1-ld) of MSPs^^ or LMXBs in 
globulars. Total DNS numbers were derived from these cross sections, NS densities and 
number of progenitor binary systems for clusters near core collapse (see text). 
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Figure 2. Binary period vs. eccentricity for DNS systems formed by NS ex¬ 
change into LMXBs with O.4M0 secondaries. Only DNS binaries produced with 
merger times ^10 Gy are plotted. Each “heavy” point marks a binary that has sur¬ 
vived until today but will merge within a Hubble time, whereas the initially created 
DNS systems are marked by the small dots (most of which have merged). The Star 
marks the parameters for the M15-C system, showing that it indeed could be produced 
by an exchange encounter of a held NS with an LMXB (or qLMXB) containing a NS and 
O.4M0 secondary. 




